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Abstract 

The absorption spectrum of the title compound in the spectral range of 
the Hydrogen-bonded OH-stretching vibration has been investigated using a 
five-dimensional gas phase model as well as a QM/MM classical molecular 
dynamics simulation in solution. The gas phase model predicts a Fermi- 
resonance between the OH-stretching fundamental and the first OH-bending 
overtone transition with considerable oscillator strength redistribution. The 
anharmonic coupling to a low-frequency vibration of the Hydrogen bond lead- 
ing to a vibrational progression is studied within a diabatic potential energy 
curve model. The condensed phase simulation of the dipole-dipole correlation 
function results in a broad band in the 3000 cm -1 region in good agreement 
with experimental data. Further, weaker absorption features around 2600 
cm -1 have been identified as being due to motion of the Hydrogen within the 
Hydrogen bond. 
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1. Introduction 



Infrared (IR) spectroscopic studies of the dynamics of Hydrogen bonds 
(HBs) continue to trigger considerable theoretical and experimental efforts 
not at least due to the outstanding importance of this type of bonding for 
a variety of phenomena in chemical and biological physics [H, 0]. In recent 
years ultrafast IR spectroscopy has demonstrated its capability to unravel 
the details of the dynamics hidden in linear absorption spectra [3|. Phthalic 
acid monomethylester (PMME) has served as model system for studying the 
ultrafast vibrational dynamics of medium strong HBs jl, 0, 0, 0, 0, 0, E3] ■ 
Experimental results include the observation of wave packet dynamics due 
to excitation of a low-frequency mode which modulates the HB length and 
which is coupled to the OH-stretching fundamental transition 0,0], and the 
unraveling of the path for the ultrafast vibrational energy relaxation using 
two-color pump-probe spectroscopy The observed subpicosecond relax- 
ation of the OH-stretching mode has been explained by an energy cascading 
mechanism involving different OH-bending modes of the HB. Concerning the 
linear IR spectrum, however, there is only a study for the deuterated case 
which combined a three-dimensional model Hamiltonian with rates for phase 



and energy relaxation according to dissipation theory [10J. Specifically, it 
was found that the IR lineshape in the OD-stretching region is dominated 
by a Fermi-resonance between the OD-stretching fundamental and the OD- 
bending overtone. The resulting double peak sha pe of the absorption is 



typical for Hydrogen-bonded systems, see Refs. [111. |12[|. 

The present study sets the focus on the IR spectrum of the normal species 
of PMME in the range of the OH-stretching fundamental transition around 
3000 cm -1 . Two complementary models will be considered. First, a five- 
dimensional gas phase Hamiltonian is determined and diagonalized to obtain 
detailed insight into the composition of vibrational eigenstates in the consid- 
ered spectral range. Degrees of freedom are selected starting from the OH- 
stretching and bending modes which are supplemented by additional modes 
according to calculated force constants up to fourth order. Second, hybrid 
quantum mechanics/molecular mechanics (QM/MM) simulations of classical 
trajectories of PMME in CCI4 solution are performed. Besides an analysis 
of the HB geometry in solution, the dipole-dipole autocorrelation function is 
calculated from which the IR absorption spectrum is obtained (for related 



applications to Hydrogen-bonded systems, see Refs. [13j, [lj, |l5|, lly, Il7l|). 



The paper is organized as follows: Section [2] starts with an account on the 
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determination of the Hamiltonian for PMME in gas phase. This is followed by 
a summary of the QM/MM protocol. Results for the gas phase IR spectrum 
are given in Section 13.11 and the analysis of the condensed phase simulations 
is presented in Section 13.31 The paper is summarized in Section HI 



2. Theoretical Methods 

2.1. Anharmonic Potential Energy Surface in Gas Phase 

The electronic ground state geometry of PMME in the gas phase has 
been optimized using the DFT/B3LYP level of theory with a Gaussian 6- 
31+G(d,p) basis set [l8[, see Fig. [TJ A detailed account on the gas phase 
geometry as well as on its dependence on the quantum chemical method 
can be found in Ref. Focussing on the vicinity of the most stable 

configuration we have chosen to express the gas phase Hamiltonian in terms 
of normal mode coordinates!^}. In general one can write the PES in terms 



of a correlation expansion [20 



V({Qi}) = X)V 4 (1) (O0+X)^ a) (Qi,Oi)+ E vfitiQitQitQk) 

i i<j i<j<k 

+ E *®Q„ Q fc , Q,) + . . . C 1 ) 



i<j<k<l 



Restricting to four-mode correlations this type of PES is readily obtained by 
the anharmonic expansion 

v({Qi}) = ^ E + j\ E K »kQiQjQk + \ { E K mQiQjQkQi ■ (2) 

i ijk ijkl 

Here, the u>i are the harmonic frequencies and the and Kijki are the third- 
and fourth-order anharmonic couplings. For the present case of PMME they 
have been calculated by combining finite differencing with analytical second 
derivatives as described in Refs. 21, 22 



Being interested in the PMME HB dynamics we have selected a set of rele- 
vant modes comprising a reduced five-dimensional (5D) quantum mechanical 
model as follows: Starting point are the coupled OH-stretching (Voh) and 
OH-bending (Son) modes with harmonic frequencies of 3279 cm' 1 and 1446 
cm -1 , respectively, shown in Fig. [2j Further we included the strongest 
coupled low-frequency mode i^hb with harmonic frequency of 39 cm -1 . As 
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can be seen from Fig. [2] excitation of this mode leads to a periodic modu- 
lation of the 0-0 distance in the HB. Hence the choice of such a mode is 
motivated by the experimental observation of damped low-frequency wave 
packet motion in Ref. Note, however, that from the oscillatory com- 
ponent of the pump-probe signal a frequency of ~100 cm" 1 had been ex- 
tracted. Looking for higher frequency modes of similar character one finds 
a vibration with a harmonic frequency of 72 cm -1 [3]. However, its anhar- 
monic coupling to the stretching vibration is considerably smaller than for 
the selected z^hb- Further, for such low-frequency modes the harmonic ap- 
proximation performs usually rather poor. Indeed the anharmonic frequency 
of z/hb within the present model is 60 cm -1 . The selection of further modes 
is guided by the available anharmonic coupling constants as well as the possi- 
ble resonances. Here, we find two modes, 71 and 72, which have out-of-plane 
OH-bending character (see, Fig. [2]) that could have significant influence on 
the OH-stretching and in particular OH-bending vibration. Their harmonic 
frequencies are 682 cm -1 and 785 cm -1 , respectively, and their combination 
as well as the overtone transitions are close in resonance to the fundamental 
transition of mode 5qh- Moreover, the combination of 5oh and the overtones 
of 71 and 72 is close in energy to the OH-stretching fundamental transition. 
The most important anharmonic coupling constants for these five modes 
are summarized in Tab. HJ The 5D model coordinates will be labeled as 

{Qs, Qb, Q71, Qhb}- 

The question arises whether the fourth-order approximation, Eq. ([2]), is 
justified for describing the selected modes. In order to scrutinize this point 
we plot in Fig. [3] exemplarily the potentials along Q s , Q^, and Qhb as ob- 
tained from the fourth-order expansion (dashed line) and from a numerical 
evaluation on a grid (solid line). According to panel (a), the anharmonic 
force field description of z/qh reproduces the exact potential for small dis- 
placements, but the two curves start to deviate above ~4000 cm -1 . The 
anharmonicity of the bending mode, 5on, is not very pronounced what is 
reflected in a good agreement between the two curves up to 12000 cm -1 in 
panel (b). The fourth order description breaks down completely for the low- 
frequency mode, z/hb, as seen i n panel (c). The anharmonicity of 71 and 
72 is only modest and fourth-order expansions and grid based potentials es- 
sentially agree over the considered range. Nevertheless, we have chosen to 
use the obtained grid-based data for representing all one mode potentials 
Vi (Qi). Further, the following two-mode potentials have been generated 
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on a numerical grid: V S ^(Q S , Q h ), V S ^ B (Q B , Q H b), and V$ B (Q h , Que)- This 
selection is motivated by the focus which is on the OH-stretching region 
of the spectrum where (i) the Fermi-resonance type coupling between the 
stretching fundamental and the bending overtone is of particular importance 
(see Tab. [[]) and (ii) low-frequency wave packet motion can be triggered by 
OH stretching excitation. All other correlations up to the four-mode term 
have been included in terms of the calculated force constants. 

The resulting Hamiltonian has been diagonalized assuming a diagonal 
kinetic energy operator in a step- wise procedure as follows: First, the un- 
coupled one- dimensional Hamiltonians have been diagonalized (note that the 
coordinates are mass-weighted) 



h2 d 2 . 



2 dQ\ 



Wn) = KWn) (3) 



by using the Fourier- Grid-Hamiltonian (FGH) method [23J (Grid parameters: 
Q s - 128 points in [-0.9:1.7], Q h - 64 points in [-2.0:2.0], Q 7l and Q l2 - 64 points 
in [-2.1 2.1], Qhb - 128 points in [-8.0:8.0] (intervals in a oy /&.m.u.)). Second, 
from these zero-order functions {|0^)}, two product bases are formed for the 
modes (t'oH, #oh) and (71, 72) including 5x7 and 7x 7 functions, respectively. 
The respective 2-mode Hamiltonians are diagonalized and the lowest 5 and 
20 eigenf unctions are used to span a product basis for the diagonalization of 
the 4-mode Hamiltonian at the equilibrium geometry of the low-frequency 
mode i/hb- Notice that the parameters for the successive diagonalization 
have been chosen such as to obtain converged eigenstates for transitions up 
to about 3600 cm -1 . 

Due to the time scale separation between the four fast modes z/oh, <5oh, 
71, and 72 and the mode ^hb, the resulting eigenvalues, which span a di- 
abatic basis, can be used for the characterization of the nature of the fast 
modes' states. Denoting the diabatic basis as {|a diab )}, these states obey the 
following Schrodinger equation: 



[ E 



[-y ^2 +^ (c) (Qfast, Qhb = 0)] |« diab ) = £f*V iab > 



i=s,b,7i,72 



2 dQi 



(4) 

Here, Qf as t comprises the four fast modes, whereas V^(Qfastj Qhb = 0) 
stands for the couplings between them. Expressed in the basis of the uncou- 
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pled single mode states the four-dimensional diabatic states read 



|a diab ) 



^<% u \fi)\fimw) ( 5 ) 

ijkl 



with the expansion coefficients C^ kl . In the simulations below we have in- 
cluded the lowest 25 diabatic states; for an assignment see Tab. [2J 

Having defined the diabatic states, the total Hamiltonian in the diabatic 
representation obtains the following form 

#diab = ^2 \^ a P ( _ Y~dCP ^ E " ab + ^(^hb)) 

+ (l-^)^(Q HB )]|« diab )(/3 diab | (6) 
where the matrix elements between the diabatic states are given by 

V^CQhb) = (« diab |V f S 1 B ) (Q H B)+^ (c) (Qfast,QHB)-^ (c) (Qfas t ,QHB = 0)|/? diab > . 

(7) 

In a final step the Hamiltonian, Eq. ([6]), is diagonalized after introducing 
zero order states of the low-frequency mode with respect to the different 
diabatic states' potential energy curves. Labeling the resulting 5D eigenstates 



as the IR stick spectrum can be expressed as 24 



/(^h oc u E p <*i d Si 2 ^ - ^) ■ ( 8 ) 

i=x,y,z a{3 

Here, P a is the thermal distribution function, and d^p and oo al 3 stand for 
the zth transition dipole moment component and the transition frequency, 
respectively. For the dipole moment vector d we have incorporated its three 
components, S l \ by expressing them on one-dimensional grids along Q s and 
Qb- 111 other words, only selected one-mode contributions to the dipole mo- 
ment function are considered. 

2.2. QM/MM Approach 

PMME solvated in CC1 4 is treated on the basis of the hybrid QM/MM 



method provided by the CPMD/Gromos interface [25j. The QM/MM sep- 
aration has been done in the solute/solvent fashion. The solvent molecules 
are coped with Gromos adopting the Gromos96 force field [26j| . The PMME 
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molecule is dealt with CPMD [27j using the Becke exchange and Lee- Yang- 



Par correlation functional (BLYP) [28|, [29|. Further, the Troullier- Mart ins 
(TM) pseudopotential [30] with a wavefunction cutoff of 70 Ry is adopted 
to describe the interaction between the core and valence electrons in the 
quantum region. For the fictitious electron mass we use 400 a.u. 

In the simulation, one PMME is solvated in 769 CCI4 molecules and 
the solution is put into a box with dimensions 50Ax50Ax50A. Before the 
QM/MM simulation, a classical solvent equilibration with Gromos is carried 
out for 1 ns at 300 K with the solute fixed using the SHAKE algorithm. The 
quantum part is placed in a 15.88 A x 12.70 A x 19.05 A box. The trajectory 
run is performed at 300 K with a Nose- Hoover chain thermostat and by using 
a time step of 2 a.u. (0.048 fs). The total simulation time is 16 ps, which 
takes about 40 days on a cluster with 24 2.66 GHz CPUs. The first 0.5 ps 
has been assigned for equilibration of the QM/MM box. 

The resulting trajectory has been analyzed putting emphasis on the ge- 
ometry of the HB. Further, dipole and velocity autocorrelation functions 
were calculated to give the IR absorption spectrum and the density of states, 
respectively. The IR absorption coefficient of the solute, / d-d (u;), can be 



calculated from the dipole-dipole correlation function 24 



"OO 

I d -'\u) oc c^Re / dte iuJt (d(t) ■ d(0)) , (9) 







whereas the density of states is given by [31 



N(u)= I die^^< Vi (i)-Vi(0)). (10) 

Further, we have calculated anharmonic OH-stretching and bending frequen- 
cies along the trajectory by generating two-dimensional snap-shot potential 
energy surfaces for an otherwise fixed geometry of solute and solvent 32j, 133 



Anharmonic vibrational eigenstates have been obtained using sine basis sets 



34] and including the kinetic energy coupling between the coordinates. The 
potential is obtained by moving the H atom on a two-dimensional plane de- 
fined by the atoms C8-01-H. In the calculation 24 snap-shots are sampled 
starting at the production part of the trajectory with an interval of 2000 
steps. For each snap-shot, the Ol-H bond stretching is sampled with 7 grid 
points in the region 0.7-2.0 A. The Ol-H bending within the C8-01-H plane 
is sampled with 7 points from -45° to 45°. 
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3. Results and Discussion 



3.1. Analysis of Gas Phase IR Spectrum 

The decomposition of the diabatic states of the fast modes is given in Tab. 
|2]for the range 2200-3500 cm -1 . Most of the states correspond to almost pure 
combination and overtone transitions. Exceptions are the states |a diab = 16) 
and |a diab = 20). Both are mixed states with respect to the bending overtone 
and the stretching fundamental transitions. They are found at 2853 and 
3044 cm -1 , respectively. Interestingly, both states have contributions from 
combinations of the Son fundamental and the 71/2 overtone transitions. 

The diabatic potential energy curves V aa (QnB), Eq. (JTj), are depicted 
in Fig. HI The potential curves for the diabatic states |a diab = 16) and 
j^diab _ (marked in black) are found to be embedded in the manifold 
of curves originating from nearby states. Since there are no symmetry re- 
strictions, the low-frequency mode in principle will couple all states in this 
region of the present model. This gives rise to the stick spectrum, Eq. (jSJ), 
shown in the insert of Fig. HI It is dominated by the above mentioned Fermi- 
resonance, but also contains the signature of the progression with respect to 
the low-frequency mode on the blue side of the z/ Q h dominated peak. Judg- 
ing the intensities of te various peaks, it should be kept in mind that the 
dipole surface is described by selected one mode terms only. Note that this 
progression is the origin of the wave packet motion of the low-frequency HB 
mode which has been observed in Ref. [§]. Needless to point to the pro- 
nounced anharmonicity of i^hb which together with the many diabatic levels 
coupled in this spectral range makes the dynamical problem considerably 
more complicated than that of shifted harmonic oscillators. 

3.2. Hydrogen Bond Geometry in Solution 

The geometry of the HB along the trajectory will be characterized by 
the HB length, -Roi-03> the angle, (3 = Z(01,H, 03), the bending an- 
gle of the OH vibration, a = Z(C8,01,H), and the dihedral angle = 
Z(C8, 01, H, 03); for labeling see Fig. [TJ The change of these parameters 
along part of the production trajectory is given in Fig. [5j Note that the OH 
bond length itself does not change appreciably, its average value is 1.02 A 
and the variation is 0.04 A. In other words, this system does not support 
any Hydrogen transfer similar to the prediction of the gas phase calculations 



19j. From Fig. [5^ we obtain the average HB length of #01-03 = 2.63A and 



a RMS variation of 0.11 A. As compared with the gas phase DFT/B3LYP 
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6-31+G(d,p) geometry the HB is on average elongated by 0.07A [lij]. Large 
elongations correlate with strong deviations of the HB from linearity as seen 
in panel (c) of Fig. [51 e.g. for ~1.5 ps and ~3.2 ps. The average deviation 
from linearity is 24°. The time evolution of the OH bending angle a is given 
in Fig. [5b, its average and RMS variation are 111° and 6°, respectively. Fi- 
nally, we focus on the planarity of the HB with respect to the dihedral angle 
<p. Note that the gas phase structure of PMME has been calculated to be 



overall nonplanar [19(. The gas phase DFT value is obtained as = —51°. In 
condensed phase the average and RMS variation are 15° and 65°, respectively. 
The distribution of <fi values is, however, bimodal with the main peak at 57° 
and a 40% smaller peak at -53°. Thus, in comparison with the gas phase 
PMME assumes more often a structure where carboxyl and ester groups are 
twisted with respect to each other. This fact also indicates the limitations 
of the gas phase model presented above. Finally, we point out that the fre- 
quently observed changes from 180 to -180° and vice versa are indicative of 
linear HB geometries. 

3.3. IR Spectrum in Solution 

The total IR spectrum according to Eq. (Q is shown in Fig. [61 Our 
focus will be on the OH-stretch region around 3000 cm -1 which is shown 
enlarged in Fig. Panel (a) compares the calculated IR spectrum with the 
experimental result of Ref. [8j. The experimental spectrum exhibits a broad 
band covering the range from ~2750 to ~3250 cm -1 . Comparison with the 
spectrum of the deuterated species [3] shows that the sharp feature around 
2950 cm -1 must be due to some CH stretching vibration which is not coupled 
to the HB motion. The two features to the left and right of this peak survive 
deuteration as well. The two peaks around ~2600 cm -1 , on the other hand, 
disappear upon deuteration and should be related to the HB. 

The calculated IR spectrum, although not reproducing the full width of 
the measured band, compares rather well with experiment. This concerns the 
position of the main band, but also the appearance of some intensity around 
2600 cm -1 . The assignment of the spectrum can be guided by comparison 
with the stick spectrum of the gas phase model which is also shown in Fig. 
[7k. Essentially, the two Fermi-resonance lines are in the relevant spectral 
region of the main experimental band, giving evidence that the latter is 
indeed containing these transitions. Its width is likely to originate from the 
fluctuations of the HB parameters as given in Fig. [5j In addition the position 
of the stick spectrum agrees nicely with that coming from the condensed 
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phase simulation. As far as the features around 2600 cm -1 are concerned we 
can only conclude that the respective transitions are not due to excitations 
of those modes which are part of the 5D model, see Fig. [2j 

Another means for gaining insight into the nature of the motions con- 
tributing to the IR spectrum is an analysis in terms of the density of states 
due to certain groups of atoms. Respective results are given in Fig. [7b. 
Overall the width of IR spectrum in the considered range is comparable to 
the density of states. The latter is, however, dominated by contributions of 
non-H-bonded Hydrogen motions. Contributions from the Hydrogen of the 
HB are in the range between ~ 2850 and 3000 cm -1 which is in accord with 
the discussion of panel (a). Interestingly there is some intensity of N(u) 
around 2600 cm -1 , i.e. the mentioned features in the IR spectrum are in- 
deed likely to be due to the dynamics of the HB as suggested by deuteration 
studies. 

The quantum mechanical calculation of the linewidth in condensed phase 
requires to have at hand the fluctuations of the transition frequencies of the 
high frequency modes due to their interaction with the surrounding modes 



351 ] . In Ref. [36J we have used an empirical mapping between HB distances 



and transition frequencies to predict the linewidth of the NH-stretching tran- 
sition in a solvated adenine-uracil base pair. This mapping was based on 



available crystal structure data [37|, see also Ref. [38|. However, no such 
empirical correlation is available for PMME type of molecules. An alterna- 
tive could be the explicit calculation of potential energy curves for selected 



coordinates on-the-fly along the trajectory 32j, |33j. Figure [8] shows results 



of respective calculations for the Fermi-resonance states, i.e. the dependence 
on the bending overtone and the stretching fundamental transition on the 
HB length, i? i-03- F° r the i^oh mode one can recognize the expected cor- 
relation, i.e. the increase of the transition frequency with i?oi-03- Note 
that this holds irrespective of the nonlinearity and wide variability of the HB 
parameters seen in Fig. [5j For the bending overtone, on the other hand, 
no such correlation is discernible. The reason can be found in the difficulty 
to define the coordinate for this bending motion due to the large flexibil- 
ity of the structure which mixes different types of bending motions, a fact 
that cannot be captured by the present definition of coordinates within the 
two-dimensional model. 
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4. Summary 



The IR absorption spectrum of PMME in the spectral range of the Hydro- 
gen-bonded OH-stretching vibration has been investigated using a five-dimen- 
sional gas phase quantum mechanical model as well as condensed phase clas- 
sical molecular dynamics simulations on the basis of the QM/MM approach. 
The gas phase model predicts specific anharmonic coupling patterns between 
the high-frequency OH stretching vibration and lower-frequency bending 
modes. In particular it was found that the region around 3000 cm -1 contains 
a Fermi resonance between the z/QH-stretching fundamental and the first <5oh- 
bending overtone vibration with considerable oscillator strength redistribu- 
tion. The effect of low-frequency HB motion has been studied within the 
picture of potential energy curves defined for the diabatic states of the higher- 
frequency modes. These curves are markedly anharmonic and form a dense 
manifold of coupled states in the 3000 cm -1 range. In the IR spectrum the 
low-frequency mode gives rise to a vibrational progression along with the ^oh 
fundamental transition. These finding are similar to the case of deuterated 
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The condensed phase simulations of the IR spectrum gave a broad band 
in the 3000 cm -1 region whose position is not only in agreement with the 
main experimental band. In addition there is a weak band around 2600 
cm -1 coinciding with an experimental feature. Analysis in terms of the den- 
sity of states led to the conclusion that this band has contributions from 
HB motions. Comparing gas and condensed phase IR spectra we found an 
overall good agreement as far as the spectral range is concerned. However, 
the trajectory simulations also indicate that in solution a different config- 
uration may also be present and contribute to the spectrum. Finally, we 
have explored the possibility to correlate i^oh fundamental and 5oh overtone 
excitation to the length of the HB. While there appears to be a reasonable 
correlation with respect to von, the scatter in the <5oh points is considerable, 
thus pointing to the difficulty of defining the bending motion in terms of a 
single coordinate. 
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Table 1: Selected cubic (left) and quartic 
chosen 5D gas phase model. 
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right) anharmonic coupling constants for the 
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Table 2: Energies of the diabatic states in the interval from 2200 to 3500 cm -1 assignment, 
expressed in terms of the uncoupled anharmonic modes. The last column contains the 
expansion coefficients (> 0.1 for states 10-21 and only leading coefficients for states 22- 
24). 
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Figure 1: Structure of phthalic acid monomethylester (PMME) as obtained by geometry 
optimization using the DFT/B3LYP level of theory with a Gaussian 6-31+G(d,p) basis 
set. In the analysis of the trajectory the angles a — Z(C8, 01, H), (3 = Z(01, H, 03), and 
the dihedral angle cf> formed by the atoms C8, 01, H, and 03 will be used. 
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Figure 3: Potential energy curves along modes i/oh (a), 5oh (b), and i^hb (c). The solid 
lines correspond to the potentials obtained on a grid, whereas the dashed lines stand for 
the Taylor expansions up to the 4 th order according to Eq. (JSJ. 
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Q HB (a (a.m.u.) 1/2 ) 

Figure 4: Diabatic potential energy curves of the lowest 20 states. The thick solid lines 
mark the states which are dominated by the bending fundamental, |a dlab = 5), its first 
overtone |a dlab = 16) as well as the stretching fundamental, |a dlab = 20), transition. The 
insert shows the stick spectrum according to Eq. {8j). 
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Figure 5: Geometric parameters of the HB in PMME along a representative part of the 
trajectory : (a) HB length, i?oi-03j (b) bending angle a = Z(C8,01,H), (c) out-of- 
line motion of the hydrogen measured by the difference between 180° and the HB angle 
P = Z(01,H,03), (d) dihedral angle of the HB (ft — Z(C8, 01, H, 03). The horizontal 
line in each panel indicates the corresponding time average over the whole production 
trajectory. 
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Figure 6: Total IR spectrum of PMME in CC1 4 solution at 300K as obtained from Eq. © 
using the 15.5 ps production QM/MM trajectory. Note that the contribution of solvent 
absorption is not included since the dipole autocorrelation function is calculated for the 
QM part only. 
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Figure 7: (a) Experimental (dashed) QM/MM (solid), and gas phase (sticks) IR 
spectra, (b) Comparison between the density of states, Eq. (fTUl) . and the IR spectrum as 
indicated in the figure key. 
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Figure 8: The correlation between the HB bond length, i?oi-03: and the OH-stretching 
fundamental transition frequency (bullets) as well as the first bending overtone (hollow 
squares). For definition of coordinate see text. 
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